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ABSTRACT 

Intermediate resolution {R ~ 7000) spectroscopy is presented for 76 photometrically 
selected very low mass (0.04 < M < 0.3 Mq) candidate members of the young cluster 
around a Orionis. More than two thirds appear to be genuine cluster members on the 
basis that they exhibit Lii 6708A absorption, weak Nai 8183/8195A features and a 
radial velocity consistent with the cluster mean. Photometric selection alone therefore 
appears to be very effective in identifying cluster members in this mass range. Only 
6 objects appear to be certain non-members, however a substantial subset of 13 have 
ambiguous or contradictory indications of membership and lack Li absorption. This 
together with an observed spread in the equivalent width of the Li absorption feature in 
the cooler stars of our sample indicates there may be deficiencies in our understanding 
of the formation of this line in cool, low-gravity objects. 

Four candidate binary cluster members are identified. Consideration of sampling 
and radial velocity measurement precision leads us to conclude that either the fraction 
of very low mass stars and brown dwarfs in small separation (a < 1 au) binary systems 
is larger than in field M-dwarfs, or the distribution of separations is much less skewed 
towards large separations. This conclusion hinges critically on the correct identification 
of the small number of binary candidates, although it remains significant even when 
only the candidate members displaying Li absorption are considered. 

Broadened Ha emission, indicative of circum(sub)stellar accretion discs is found 
in 5 or 6 of the candidate cluster members, 3 of which probably have substellar masses. 
The fraction of accretors (10 ± 5 per cent) is similar to that found in stars of higher 
mass in the a Ori cluster using Ha emission as a diagnostic, but much lower than 
found for very low mass stars and brown dwarfs of younger clusters. The timescale 
for accretion rates to drop to < 10~^^ Mq yr"-'^ is hence less than the age of the a Ori 
cluster (3 to 7Myr) for most low-mass objects. 

Key words: techniques: radial velocities - techniques: spectroscopic - stars: low- 
mass, brown dwarfs - stars: pre-main-sequence - open clusters and associations: indi- 
vidual: a Orionis. 



1 INTRODUCTION 

The typical Jeans mass in molecular cloud cores is an order 
of magnitude more massive than a brown dwarf, yet some- 
how these substellar runts are found in great abundance - 
either isolated in the field, as companions to more massive 
objects or as members of young clusters (e.g. Bouvier et 
al. 1998; Kirkpatrick et al. 2000; Bejar et al. 2001, Gizis 
et al. 2001; Moraux et al. 2001). Despite this, there is no 
universally accepted formation scenario which explains the 
profusion of substellar objects. Some theories propose that 
turbulent fragmentation can extend the more conventional 
star formation scenario to low masses (Padoan & Nord- 



lund 2002), whilst others have proposed that BDs can form 
through instabilities in circumstellar or circumbinary discs 
(Pickett et al. 2000; Jiang, Laughlin & Lin 2004); by early 
ejection from protostellar aggregates (Reipurth & Clarke 
2001); or by removal of accretion envelopes by collision or 
photoevaporation (Price & Podsialowski 1995; Kroupa & 
Bouvier 2003). Possible diagnostics that may point to which 
processes are the most important include the initial mass 
function (IMF) and the properties of binary systems and 
discs among the lowest mass objects. 

The IMF is often represented by the power-law relation- 
ship; dn/d\ogm oc m~". Since the classic work of Salpeter 
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(1955) the IMF has been well studied down to masses ol 
order ~ 1^0 ^^id is best described by an exponent ol 
Q = 1.35. Below this it flattens, with q ~ having been 
measured in the field and in several young open clusters 
(see Chabrier 2003 for a review). For VLMS and BDs the 
IMF is more uncertain. Most determinations in young clus- 
ters report —0.5 < a < (e.g. Moraux et al. 2001; Bejar 
et al. 2001), which are consistent with determinations in 
the field once appropriate corrections for unresolved bina- 
rity are applied (Chabrier 2003). Recently there have been 
claims for non-universality of the IMF at the lowest masses. 
Briceno et al. (2002) and Luhman et al. (2003) report a sig- 
nificantly smaller number of BDs per star in Taurus- Auriga 
and IC 348 than in the Orion Nebula cluster, perhaps indi- 
cating an important environmental dependence for the BD 
formation process. 

Determination of the true IMF is inextricably linked 
with a knowledge of binary properties. VLMS and BDs can 
be hidden in unresolved binary systems, changing the ob- 
served mass functions and possibly introducing apparent 
IMF variations between clusters with different binary fre- 
quencies (Kroupa & Bouvier 2003). The statistics of binary 
systems themselves - mass ratios (g) and separations (a) ~ 
are also important. The binary properties of BDs and VLMS 
are not simply an extension of those for low-mass stars. The 
lack of wide binary VLMS and BD systems (a > 20 au) in 
the field and young clusters is notable (Martm et al 2003; 
Bony et al. 2003; Close et al. 2003; Gizis et al. 2003) and has 
been taken as support for dynamical ejection models. The 
overall binary frequency is still uncertain. Directly resolved 
binary systems with 1.6 < a < 16 au and q ^ 0.5 form about 
15 per cent of the population. This figure is reasonably con- 
sistent with, although a little higher than, the number of BD 
binaries with separations > 1 au formed in hydrodynamical 
simulations (Bate, Bonnell & Bromm 2003). Closer binaries 
are much less well studied. Examples of spectroscopic BD 
binaries exist (e.g. PPL 15 - Basri & Martin 1999) and may 
be quite common (Reid et al. 2002; Guenther & Wuchterl 
2003). An indirect estimate of the fraction of unresolved bi- 
naries in the Pleiades, using position in colour-magnitude 
diagrams, led Pinfield et al (2003) to conclude that the BD 
binary fraction with q > 0.4 in the Pleiades was 50^itJ per 
cent. The implication is that most of these must be close 
binaries (a < 1.6 au) in order to agree with the statistics 
for wider binary systems found in high spatial resolution 
imaging. 

Insights into VLMS and BD formation may also come 
from the presence (or otherwise) of circumstellar discs, 
which can be characterised by strong Ha emission or in- 
frared excesses. Bate et al. (2003) argue that the ejection 
of a brown dwarf will result in the truncation of any sur- 
rounding material leading to less massive discs with shorter 
lifetimes. However, recent results using accretion diagnos- 
tics and I/'-band excesses suggest that many BDs do indeed 
possess (accretion) discs and that their lifetimes may be sim- 
ilar to discs around more massive objects (see Liu, Najita 
& Tokunaga 2003; Jayawardhana, Mohanty & Basri 2003; 
Jayawardhana et al. 2003b). 

To a greater or lesser extent, studies of the IMF, bina- 
rity and discs in young clusters are hampered by a lack of 
firm membership information, and in the case of close bi- 
narity in VLMS/BDs, a lack of data. Samples must be free 



of contamination to draw accurate conclusions about the 
IMF and disc frequencies. Here we explore these issues using 
the VLMS/BD population clustered around the 09.5V star 
a Orionis, which is itself a member of the Orion OB lb asso- 
ciation (Walter et al. 1998). At an Hipparcos distance of only 
352t^f pc, aged 3 - 7Myr (Oliveira et al. 2002) and with 
low reddening {E{B — V) = 0.05 mag. Brown et al. 1994), 
the a Ori cluster is a favourable environment in which to 
study VLMS/BD evolution. Its youth guarantees that low- 
mass members will be relatively luminous. The VLMS and 
BD population of a Ori has been investigated by a number 
of authors. Bejar et al. (1999, 2001) have combined optical 
and near-infrared photometric surveys of the a Ori region to 
locate candidate VLMS and BDs and estimate an IMF with 
Q = —0.2 ± 0.4. There is some evidence that the IMF ex- 
tends down to a few Jupiter masses (e.g. Zapatero-Osorio et 
al. 2000). Low resolution spectroscopy, with resolving power 
250 < R < 500, has been obtained for many of these objects 
(Bejar et al. 1999; Zapatero-Osorio et al. 2000; Barrado y 
Navascues et al. 2003). Almost all are claimed to be mem- 
bers based on their spectral types, Ha emission character- 
istics and for some objects, on the basis that their neutral 
alkali lines (Ki, Nai) are weaker than expected for dwarf 
stars, indicating a low gravity. In our opinion, only the latter 
can reasonably discriminate against contaminating objects 
in photometrically selected samples. A spectral type consis- 
tent with colour is not evidence for membership because the 
dominant contaminants will be field M dwarfs with similar 
spectral types. Nor is moderate Ha emission conclusive. A 
large fraction of field M-dwarfs show Ha emission, especially 
in the range M3-M6V (Gizis, Reid & Hawley 2002). 

Higher resolution spectroscopy offers a number of ad- 
vantages in determining membership credentials. The pres- 
ence of Li in the atmosphere can be diagnosed from the 
Li I 6708A line, indicating a youthful star. This technique 
works because Li is burned very quickly in fully convective 
VLMS envelopes once the core reaches ~ 3 x lO" K. Further- 
more, radial velocities (RVs) allow the kinematic properties 
of candidates to be used as a check on membership. Lit- 
tle work at high resolution has been done on VLMS and 
BDs around a Ori. Zapatero-Osorio et al. (2002) obtained 
spectroscopy of 26 low-mass stars and 2 BD candidates at 
1600 < R < 4400, of which only 5 had / > 15, corresponding 
to masses < O.2M0 . These data demonstrated that Li should 
remained unburned in all VLMS and BD cluster members. 
MuzeroUe et al. (2003) obtained R = 8000 and R = 34000 
spectroscopy of 7 BD candidates from Bejar et al. (1999, 
2001) with 16.5 < / < 19.8, finding that 2 have radial veloc- 
ities inconsistent with cluster membership and showing no 
Li. Finally, McGovern et al. (2004) observed 2 photometric 
BD candidates at i? = 1000 — 2000 in the near infrared, 
showing that one object had neutral alkali lines too strong 
to be a pre-main-sequence (PMS) object. 

In this paper we present new RI photometry and R ~ 
7000 optical and near IR spectroscopy of 76 photometrically 
selected VLMS/BD candidates around a Ori with 14.8 < 
I < 18.2. Our objectives are as follows. 

(i) To confirm membership of the association for these ob- 
jects, using radial velocities and the strength of Li I and Na I 
lines, to validate IMFs derived solely from photometrically 
selected members (e.g. Bejar et al. 2001). 
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(ii) To search for short period binary systems to see 
whether there are a large number of VLMS and BDs hidden 
in such systems as suggested by Pinfield et aL (2003). 

(iii) To look for the presence of accretion betrayed by the 
presence of broad Hq emission lines and compare disc fre- 
quencies in the VLMS/BDs of a Ori with those in stars of 
higher mass and with objects of similar mass in younger 
clusters. 

Section 2 of the paper outlines the optical dataset used 
to select candidate cluster members; section 3 describes the 
multi-object fibre spectroscopy; section 4 presents the spec- 
troscopic results and radial velocities; cluster membership 
is discussed in section 5 and in section 6 we discuss the 
use of various membership indicators in star forming regions 
and the IMF, binarity and accretion statistics of the cluster 
around a Ori. Our conclusions can be found in section 7. 



2 OPTICAL PHOTOMETRY 

We observed the field around a Ori with the Wide Field 
Camera (WFC) on the 2.5-m Isaac Newton Telescope using 
Harris R and Sloan i filters on the nights of 27 — 30 Septem- 
ber 1999. The WFC consists of a mosaic of 4 2048x4096 
pixel CCDs covering a sky area of about 970 square ar- 
cminutes and we obtained 8 overlapping pointings, only 5 
of which (pointings 1, 2, 3, 4 and 8) are relevant to this 
paper. The details of our observations are given in Tabled 

Each flatfield, bias and data frame was bias subtracted 
using the mean of an over-scan region, and then linearised 
using the September/October 1999 coefficients. After this, 
a median bias frame was subtracted from the flats and data 
frames to remove any residual bias pattern. Since the cam- 
era was suffering from light leaks during the run, the twi- 
light flats could only be used to remove small scale pixel- 
to-pixel variability. The large scale flatflelding was achieved 
using flats taken in March 1999. Finally all the i band data 
were defringed using master fringe frames constructed from 
data taken in February 1999. These frames were created 
by normalising and then median stacking data frames from 
each night, and then creating an average from several nights. 
Varying fractions of the fringe frame were then subtracted 
from the data until the mean average deviation in the sky 
of the data frame was minimised. From this point the data 
reduction followed the procedure described in detail in Nay- 
lor (1998) and Naylor et al. (2002). We created an initial 
catalogue of objects by searching the sum of the long i-band 
exposures, and then carried out optimal photometry in each 
separate CCD image at those positions. 

Our astrometry was initially performed using a nine pa- 
rameter model, which included the tangent point and pin- 
cushion distortion co-efficient as free parameters. Although 
we finally used the SuperCosmos sky survey (Hambly et al. 
2001) as our reference catalogue, the stars were so sparse in 
the parts of pointing 8 (see Table 1) in the direction of the 
molecular cloud, that for this pointing only, our initial solu- 
tions were with respect to the 2MASS survey (Cutri et al. 
2003). We found the mean of the distortion co-efficients de- 
rived from each CCD image was close to the nominal value, 
so we fixed its value at the nominal 220. We then fitted 
again to find the tangent points for each CCD image, which 
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Figure 1. I,R-I colour-magnitude diagram. Spectroscopic tar- 
gets are shown as follows; circles represent objects which display 
lithium in their optical spectra, crosses are definite cluster non- 
members, and our "grey area" objects are indicated by open tri- 
angles (see section 6.1). Two isochrones are shown, aged 5 and 
120Myr, derived using the models of Baraffe et al. (1998) and 
assuming a distance of 352 pc and E{R — I) = 0.03. 



we then averaged to find a tangent point for each point- 
ing. With these parameters fixed, final astrometric solutions 
with respect to SuperCosmos (with just six free parameters) 
were found. The RMS deviation from the astrometric solu- 
tions for each CCD ranged between 0.13 and 0.20 arcsec. 

To correct our optimal fiuxes to those that would have 
been obtained from a large aperture we used second and 
third order polynomials to model the profile-correction in 
the directions of the short and long axes of the CCDs respec- 
tively (see Naylor et al. 2002) . The photometric co-efficients 
used to convert the profile corrected magnitudes to appar- 
ent magnitudes were derived from observations of over 700 
standard stars in R and i during the 4 photometric nights 
of the run (Landolt 1992). These were used to obtain ex- 
tinction co-efficients (1 per night), colour co-efficients (1 per 
CCD) and zero points (1 per CCD per night). However, the 
small range of airmasses we observed at (1.1 to 1.5) meant 
that we fixed the extinction at the La Palma mean values 
(0.04 in /, 0.05 in R — I). We found we had to add a large 
systematic error (0.04 mag in I 0.06 mag in R — I) to obtain 
a reduced of 1 for the fit. After doing this, our final num- 
ber of standard stars used was 653 (8 clipped out) in i? — / 
and 720 (5 clipped out) in /. Note that although we used 
a Sloan i filter, the magnitudes we derived are tied to the 
Landolt standards and so they are on the Cousins system. 

The resulting apparent magnitudes could then be com- 
pared for stars which fell in more than one CCD field. From 
the stars in each overlap region we calculated a mean dif- 
ference for that overlap. We then adjusted the zero points 



© 2004 RAS, MNRAS 000,0-?? 



4 M.J. Kenyan, et al. 



Table 1. Wide Field Camera observing log. Coordinates refer to the rotator centre of the telescope. 





Pointing 1 


Pointing 2 


Pointing 3 


Pointing 4 


Pointing 8 


RA(J2000) 
Dec(J2000) 


05 39 50 
-02 20 00 


05 39 50 
-02 51 40 


05 37 42 
-02 20 00 


05 37 44 
-02 51 40 


05 41 20 
-02 27 50 


Filter 


R i 


R i 


R 


i 


R 


i 


R i 


1999-09-27 
1999-09-28 

1999-09-29 
1999-10-01 


2x300s 2xl50s 
lx20s IxlOs 


2x300s 2x1 50s 
lx20s IxlOs 


lx20s 
lx300s 


IxlOs 
2x1 50s 


lx20s 
lx300s 
lx300s 


IxlOs 
2x1 50s 
lxl50s 


lx20s IxlOs 



3x300s 3x1 50s 



to minimise the differences, obtaining a root-mean-square 
(RMS) discrepancy of 0.007 and 0.008 mags in / and i? - / 
respectively. Where a CCD had no overlaps, its photome- 
try was adjusted by the mean correction for that pointing, 
to maintain consistency. Finally, the mean magnitude was 
calculated for stars appearing in two or more CCD fields. 

As explained in Naylor et al. (2002), the overlaps are a 
powerful diagnostic of the accuracy of the photometry, since 
they test errors in profile correction and large-scale fiatfield- 
ing, in addition to any change in transparency. They rep- 
resent a very good measure of the likely photometric error. 
In this case, before we adjusted the zero points, we found 
a RMS difference of 0.05 mags between the overlaps. Since 
this closely matches the systematic error we had to add to 
the standard stars, we believe this represents our photomet- 
ric accuracy. We regard this as rather poor, and are unsure 
of why this is so. Some of the error is due to differences be- 
tween our filters and those used to create the standard star 
system. Specifically, we have used a Sloan i filter, whilst the 
Landolt standards are for /. However, this is unlikely to be 
the entire cause of the problem, since a linear fit of _R — J vs 
R - iiov the stars of Smith et al. (2002) gives an RMS of 
0.026 mags. In other work, using the same software, we have 
found much smaller differences between the magnitudes of 
stars measured in different CCD pointings (Naylor et al. 
2002; Burningham et al. 2003). This suggests to us that the 
remaining problem originates in the flatfield, although we 
cannot prove that it is not the profile correction which is 
at fault. Our defringing technique did not entirely remove 
the fringes, leaving structure on a similar scale to our sky 
background boxes. This means the sky subtraction is not 
perfect, which increases the uncertainties in our photometry 
beyond those given by the optimal extraction. Experiments 
with photometry of empty regions of sky showed that the 
uncertainties we quote should be increased by a factor of 
about 1.5. 

We compared our catalogue with that of Bejar et al. 
(2001). After applying a shift of 1.3 arcsec in RA and 1.1 
arcsec in declination to our catalogue, we found good corre- 
lations out to 5 arcsec, with an RMS separation of 1.8 arc- 
sec between the two catalogues. Our small astrometric RMS 
with respect to SuperCosmos suggests these numbers repre- 
sent systematic shifts and uncertainties in Bejar et al.'s as- 
trometry. Choosing just the objects with a signal-to-noise ra- 
tio greater than 10 in both catalogues, we found I — I Bejar ~ 
-0.12 ± 0.04 and {R - I) - {R - I) Bejar = -0.07 ± 0.04. 



There were no significant colour terms, but the range of 
colours used is small (just the PMS). Of more concern is 
the scatter in the relationships, with the RMS difference be- 
tween measurement's of the same star in the two catalogues 
being 0.24 magnitudes in / and 0.13 magnitudes in R — I. 
Of course some of this scatter is likely due to variability as 
most of the comparison stars are probably PMS objects that 
are young, rapidly rotating and magnetically active, possibly 
with accretion discs (see Scholz & Eisloffel 2004). 



3 SPECTROSCOPY 

3.1 Target selection 

Our spectroscopic targets have been selected on the basis 
of their position on the I,R — I colour-magnitude diagram 
(CMD), see Figure We chose objects that were reason- 
ably close to a 5Myr isochrone, generated from Baraffe et 
al. (1998) models, using an empirical colour-effective tem- 
perature relationship derived from Pleiades data and an as- 
sumed distance of 352 pc (see Jeffries, Thurston & Hambly 
2001 for details of the technique). The corresponding objects 
are highlighted on a Digitized Sky Survey (DSS) image in 
Figure 121 Targets with 14.8 < I < 18.2 (corresponding ap- 
proximately to masses 0.04 < M < 0.3 Mq using a fixed 
age of 5 Myr, a distance of 352 pc and the Baraffe et al. 
1998 models) were chosen from the catalogue with the aim 
of maximising the number of targets close to the cluster 
isochrone that could be observed in a pair of multi-fibre se- 
tups (see below). Note that because of the limited number 
of spectrograph fibres and other geometric restrictions (see 
below) we were unable to place fibres on all the good cluster 
candidates in the field of view (see also section 6.3). 

Unfortunately IR photometry was unavailable at the 
time at which the spectroscopy was performed, but was later 
obtained for all our targets from the 2MASS all-sky point 
source catalogue (Cutri et al. 2003). 

3.2 Spectroscopic Observations 

Spectra were obtained with the Wide Field Fiber Optic 
Spectrograph (WYFFOS) mounted at the Nasmyth focus 
of the 4.2-m William Herschel Telescope during the nights 
of 11 and 12 December 1999. An echelle grating and order 
sorting filters were used to obtain spectra for targets in a 
field of diameter just less than 1 degree (see Fig.l^J, covering 
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Figure 2. Digital Sky Survey image of a Orionis. Tlie positions of spectroscopic targets are marked. Tlie symbols have the meanings 
described in the caption to Figure 1. 



wavelength ranges of 6390-6810A and 7820-8460A and at 
dispersions of 0.43A pixel" ^ and 0.63A pixel"'^ respectively. 
Measurements of arc-line spectra showed that FWHM res- 
olutions of ~ 0.9 and 1.3 A were attained with these two 
settings. WYFFOS had about 100 active fibres which could 
be allocated, but restrictions on the fibre placement and fibre 
proximity meant that not all of these could be used. Differ- 
ent fibre/target configurations were used on the two nights, 
with 19 (predominantly faint) objects common to both con- 
figurations. The optical fibres had a diameter of 2.7arcsec 
and to estimate the co-temporal sky spectrum we allocated 
15 fibres in each configuration to blank sky regions. 

A log of the spectroscopic observations is given in Ta- 
ble |5| In addition to several hours of target exposures on 
each night we obtained observations of a tungsten lamp to 
aid in fibre location and fiat fielding, Copper plus Neon arc 
lamp spectra for wavelength calibration at the end of each 
target exposure, several offset exposures of blank sky that 
were median stacked to calibrate the fibre-to-fibre transmis- 
sion efficiencies in each configuration and wavelength set- 
ting, spectra of a number of M dwarfs that can be used as 
spectral-type and radial velocity standards and a spectrum 
of a bright, B3V star (HR 5191) that was used to identify 
and correct for telluric absorption lines. 



3.3 Data Reduction 

Images were processed using several packages from within 
the IRAF^ environment. Most notable of these was wyf_red 
(Lewis 1996), a script based on the fibre reduction pack- 
age DOFIBERS, but with some differences. In particular, 
WYF_RED performs CCD processing itself if provided with 
bias frames and flat fields. Complete reduction involved sev- 
eral additional tasks, the most important of which are de- 
tailed below. 

(i) Cosmic rays were removed automatically from target 
frames but a manual inspection of the arc lamp exposures 
was required to avoid the algorithm mistaking arc lines for 
cosmic rays. 

(ii) Scattered Hght was subtracted by fitting a smooth 
function to the Hght found between the fibres and then in- 
terpolating this to the regions occupied by the fibre spectra. 

(iii) The spectra were extracted using an optimal algo- 
rithm which gave more weight to regions of the profile with 
relatively high signal-to-noise ratio (SNR). 

^ IRAF is distributed by National Optical Astronomy Observa- 
tory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation. 
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6690 6700 6710 6720 6730 6740 
Wavelength (A) 

(a) Sample of optical spectra displaying the Li I line at 6708A. 
Spectra are normalised within the region 6700 — 6715A prior 
to a constant offset being applied, labels are as given in Ta- 
ble |H] The gaps in the spectra correspond to a poorly sub- 
tracted bright Sll sky line at 6717A. 
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(b) Near-infrared spectra showing the 8183,8195A Nal dou- 
blet. The dashed spectra are of the M6V field dwarf GL412b 
to demonstrate the difference in equivalent widths. Spectra 
are normalised using the region 8160-8220A and a constant 
offset applied, labels are as given in Table IHl 



Figure 3. Sample optical and near-infrared spectra. 



(iv) Each fibre has a different transmission efficiency. 6690 — 6714A region tlie SNR (per pixel) ranges from 3 to 
WYF_RED uses the offset blanlc slcy spectra to normalise tiie 20, with a median of 8. Similarly spectra in the 8175 — 8210A 
fibre responses and also correct for vignetting. window had a SNR range of 4 to 21 with a median of 9. 

(v) Copper-neon arc lamp spectra were used to determine 
a (fourth order) polynomial dispersion function for each fi- 
bre. The RMS residuals to the arc line fits were typically 
found to be 0.005A and 0.02A in the optical and near- 
infrared regimes respectively. 

(vi) All the derived spectra for a target obtained on one 
night were co-added. If an object was observed on both 

nights then these spectra were also co-added to find the One anomaly which requires comment is that during 

equivalent widths (EWs) of diagnostic lines, but the sepa- the reduction we found that two of the fibres allocated to 

rate spectra from both nights were retained to obtain two blank sky were seen to yield spectra with strong lithium ab- 

independent radial velocities (see section EJ. sorption features! We surmised that there was an error in 

the translation table which relates allocated fibre numbers 

A sample of the reduced optical and near-infrared spec- to spectrograph apertures. This hypothesis was confirmed 

tra are given in Figs. |3(a)1 and |3(b)| respectively. SNR per by examining archival WYFFOS data and shown to apply 

pixel was empirically estimated by measuring the residuals to five pairs of fibres. We have subsequently been able to 

to first order polynomial fits to pseudo-continuum windows. remove this source of error from our observations, and men- 

These values are conservative, as we know that many small tion it here for the benefit of those who intend to make use 

molecular lines are probably present in these regions. In the of AF2/ WYFFOS data (large fibres bundle) in the future. 
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Table 2. WYFFOS observing Log. Different fibre eonfigurations were employed on the two nights, indicated by the subscripts. 
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Figure 5. Heliocentric radial velocity (RV) as a function of I — J colour, cr Orionis non-members are marked through with crosses, 
"grey area" objects are triangles (see section 6.1). The 19 objects which were observed during both runs have the two velocities indicated 
(shown as dashed error bars); those for which only one velocity is available are shown as solid lines. The error bars are la values, the 
cluster mean is indicated by a horizontal line at velocity 31.2kms~^. 
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Figure 4. Errors introduced into radial velocity measurements 
by cross-correlating against a specific template spectrum. Resid- 
uals are the difference between the velocity derived using a spe- 
cific template and the mean velocity from all five templates. The 
dashed line represents a residual of zero in each panel. 



Table 3. M-dwarf spectral type standards. R—I and I — J colours 
are taken from Leggett (1992), heliocentric velocities from Nide- 
ver et al. (2002). 



Object 


Sp. 


R-I 


I-J 


RVj/ 


(GL) 


Type 






(kms^-"-) 


273 


MS. 5 


1.55 


1.49 


18.2 ± 0.1 


285 


M4.5 


1.69 


1.61 


26.5 ± 0.3 


388 


MS. 5 


1.42 


1.35 


12.4 ± 0.1 


406 


M6 


2.18 


2.33 


19.5 ± 0.1 


412b 


M6 


2.10 


1.97 




447 


M4 


1.68 


1.65 


-31.1 ± 0.1 



4 RESULTS 

4.1 Radial Velocities 

Radial velocities were calculated by cross-correlating the re- 
gion 8150-8250A , which included the prominent Nal ab- 
sorption doublet (although it is weaker than the standard 
stars in almost all our target spectra - see section 14.211 . 
against five radial velocity standards (see Table for de- 
tails) . This spectral window does contain telluric absorption 
which could potentially alter derived radial velocities and 
Nal EWs. We removed this contamination to first order by 
creating a calibration spectrum from the assumed smooth 
continuum of the B3V star we observed. The calibration 
spectrum was scaled according to the average airmass of 
observation (using the FIGARO program bsmult) and used 
to multiply the target spectrum. After this correction, very 
little telluric contamination could be seen. 

The spectra were cross-correlated using the iraf pack- 
age FXCOR (Kurtz & Mink 1998). The spectra were rebinned 



to a log-linear dispersion and then continuum subtracted. 
The ends of each spectrum are forced to zero with a co- 
sine bell function prior to a fast Fourier transform being 
performed. A filter is applied to the Fourier transforms to 
limit the frequency window used in the cross-correlation. 
The measured cross-correlation lags gave radial velocities 
which were then heliocentrically corrected according to the 
time of observation and the radial velocity of the standard 
(see Table EJ. 

To check the success of the cross-correlation process, 
we investigated the effect that using a specific standard star 
had on: 

(i) the cluster mean, 

(ii) the scatter of objects about this mean. 

This was done by cross-correlating all objects against each 
template spectrum after which, a mean velocity per object 
was calculated. The deviation from this mean when each of 
the templates was used is plotted as a function of I — J (a 
good proxy for spectral type) in Fig. 2] 

The top panel suggests that at the time of observa- 
tion, GL447 had a velocity discrepant by about 2kms~^ 
from that given in Nidever et al. (2002). It is also clear that 
there is more scatter in the RVs determined from GL 285. 
In what follows we use the average RV obtained from cross- 
correlation versus the spectra of GL273, GL388 and GL406 
as there seems to be little or no dependence of derived radial 
velocity on spectral type. 

Radial velocity errors were estimated by FXCOR. We 
tested their validity by generating three-hundred template 
spectra which were perturbed and degraded to simulate 
SNRs per pixel of 7, 15 and 25. These simulated spectra 
were cross-correlated with the original spectrum from which 
they were created. This allowed us to make an estimate of 
the error introduced by varying SNRs. From these simula- 
tions we concluded that the FXCOR uncertainty estimates 
were too small by a factor of ~ 1.75, which was approxi- 
mately independent of SNR. To remedy this, all the cross- 
correlation RV errors were increased accordingly. An esti- 
mate of any remaining error, as a result of spectral type 
mis-matches between the object and template spectra was 
determined from Fig.|l] The mean standard deviation in the 
RVs derived from each standard of 0.9kms~^ was added in 
quadrature to all RV uncertainties prior to cluster member- 
ship determination. External errors, which must be consid- 
ered if comparing our results with other published data, can 
be estimated from the deviations of the standard star RVs 
from their published values when their RVs are determined 
by cross-correlation against the other three standards (GL 
447 is clearly discrepant and excluded). We estimate that 
any external error is less than lkms~^. 

Radial velocities for all objects are listed in Table|S|and 
displayed in Figure |3 as a function of / — J colour. Those 
objects with two independent RV measurements from differ- 
ent nights are shown with dashed lines. We have observed 
one object (target 60 = SOri 27) in common with Zapatero 
Osorio et al. (2002), three (targets 39,55,75 = SOri 12, 29, 
40) in common with MuzeroUe et al. (2003), and one (tar- 
get 44) which appears in Burningham et al. (2004). Table |1| 
compares the velocities determined by these authors with 
the current work. Our results concur with the literature in 
the cases of targets 44, 55 and 60 whereas both targets 39 
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Table 4. Radial velocity comparisons between this work [K04] 
and values in Zapatero Osorio et al. (2002) [Z02] , MuzeroUe et al. 
(2003) [M03] and Burningham et al. (2004) [B04]. Objects are 
identified as in Table 6. 



ID 



Radial Velocity (kms ^) 
Z02 M03 B04 



K04 



39 
44 

55 
60 
75 



35.5 ± 10 



29.8 ± 0.7 



27.1 ± 1.' 



32.5 ± 3.3 



31.1 ± 4.1 



37 ± 2 
34 ± 3 
28 ± 4 
47 ± 5 
55 ± 9 



I 
> 




Figure 6. Lithium I (6708A) equivalent width versus I — J colour 
for our objects. Where no lithium line was observed, 2(t upper 
limits are shown. 



and 75 display evidence for binarity in the form of a variable 
radial velocity. However, for consistency, they are classed as 
single from the data in this paper (see section 5.1). 



4.2 Lithium & Sodium 

The EWs of both the Li I 6708A line and Nal near IR dou- 
blet at 8183, 8195A can be used as indicators of age and 
hence membership in samples of VLMS and BDs (see sec- 
tion |^. A low-order polynomial function was used to nor- 
malise the spectra using pseudo-continuum regions around 
both of these features (and estimate the SNR). The absorp- 
tion lines were modelled as simple Gaussian functions (two 
Gaussians in the case of the Nal doublet) and the results 
integrated to estimate the EWs. 1-sigma uncertainties in the 
EW measurements were estimated as 




Figure 7. Sodium I (8183,8195A) equivalent widths versus I — J 
colour. Filled diamonds indicate our standard star measurements, 
open diamonds represent M-dwarfs from Xu (1991) and triangles 
are M-dwarf measurements from Martin et al. (2004). Objects 
lying below and to the right of the dashed line are considered to 
have gravities lower than field M-dwarfs. 



SEW 



VFWHM X p 
SNR 



(1) 



where FWHM and p are the full width half maxima and 
pixel size respectively, in units of A. The pseudo-continuum 
regions either side of the lines are large enough that the 
statistical errors due to the uncertainty in the continuum 
level are negligible. The Nal EWs for the field M-dwarf 
(RV) standards were also measured in the same way. 

Figures |S| and [7| show the lithium and sodium EWs 
plotted as a function of I — J colour. Where no significant 
Li I feature was seen, 2-sigma upper limits to the equiva- 
lent width have been calculated based on a the median line 
FWHM and the SNR of the spectrum in question. This al- 
lows the distinction to be made between poor SNR spectra 
which may be masking the lithium line and objects which 
certainly do not possess the line. It is worth emphasizing 
that although bright forbidden lines of S ll produce a very 
noisy sky subtraction at 6716/6731 A, the sky background at 
±5 A from the Li I 6708 A line is low and well behaved (see 
Fig. 3). 

Our survey includes a handful of objects which have 
had their Li I and Nal features measured by other authors 
(Bejar et al. 1999; Zapatero-Osorio et al. 2002; MuzeroUe et 
al. 2003). A comparison of the EWs is presented in Table |^ 
There is reasonable agreement between our Li I EWs and 
those found in the high resolution spectra of MuzeroUe et al. 
(2003). However, there are some clear discrepancies between 
our EW values and those found in lower resolution spectra 
by Bejar et al. (1999) and Zapatero-Osorio et al. (2002). 
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Figure 8. Ha velocity profiles of eandidate acereting objeets. The top panel shows a typical sky profile (dotted line), plus a Gaussian 
profile (dashed line) with equivalent width 25A and a full width at 10 per cent of peak of itl35kms~^. Each subsequent panel also 
shows the observed Ha line (with continuum normalised to unity) , with an overplotted sky spectrum (scaled to the peak Ha flux) and a 
minimal accretion profile (the dashed Gaussian) to illustrate that the profile is (a) broader than any possible contribution from the sky 
and (b) broader or stronger than the minimum expected profile from any accreting object. 



Table 5. Sodium I and Lithium I equivalent widths for objects 
common to Bejar et al. (1999) [B99], Zapatero Osorio et al. (2002) 
[Z02], MuzeroUe et al. (2003) [M03], Burningham et al. (2004) 
[B04] and the current work [K04]. Objects are identified as in 
Table 6. Uncertainties for the equivalent widths are ±lA for the 
[B99] sample, ±0.09A for the [Z02] point and ±0.39A for the 
[B04] object. 







EW(N 


al ) 




EW(L 


I ) 


ID 


B99 


B04 


K04 


Z02 


M03 


K04 


39 


2.5 




2.23±0.20 




0.6 


0.54±0.06 


44 




2.65 


2.54±0.18 








55 


2.2 




1.82±0.20 




0.6 


0.32±0.14 


60 


^ 1.0 




2.86±0.20 


0.74 




0.31±0.06 


68 


^ 1.0 




3.68±0.31 








75 


< 1.0 




1.15±0.60 




0.5 


0.66±0.20 



4.3 Hq 

Strong Ha emission is a characteristic of young stars, both 
as a result of magnetic chromospheric activity and from con- 
tinuing accretion. Our attempts to study Ha in the a Ori 
candidates are hampered by the presence of strong, spatially 
varying Ha emission that can be clearly seen in our _R-band 
images and is presumably associated with hot gas in the line 
of sight towards a Ori. The signal from the background Ha 



is very strong and, in our 2.7 arcsec diameter fibres, swamps 
any Ha signal from the targets, making a measurement of 
the Hq EW impossible. However, the background emission 
is unresolved in our sky spectra, with a FWHM of about 40- 
50kms~^. White and Basri (2003) have demonstrated that 
Ha emission lines formed by accretion processes should be 
much broader than this. Classical T-Tauri stars typically ex- 
hibit Hq lines with a full- width at 10% of the peak emission 
of > 270kms"^ So, whilst we cannot measure the strength 
of the Hq line, we can certainly search for emission wings 
extending to more than ±135kms~^ from the line centre, 
which are uncontaminated by the background Ha emission. 

The Ha EW expected from accreting objects which 
demonstrate extended Hq emission with spectral types M3- 
M7 (the approximate range of our sample) is greater than 
25A and more usually several times stronger (e.g. see White 
& Basri 2003; MuzeroUe et al. 2003). Given the dispersion of 
our spectra and assuming the Ha profile is pseudo-Gaussian 
it is easy to show that such a line would produce a signal 
of strength > 0.7 times the continuum level at ±135kms~^ 
from the line centre. Such a signal should be easily detectable 
if present, even in our poorest spectra. 

Analysis of our spectra yields five objects which clearly 
exhibit extended Ha emission on this basis (see Fig. 8) and 
a further one (target 63) which has a marginally broadened 
line which may simply be a case of poor sky subtraction or 
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the result of small changes in the instrumental dispersion 
profile over the detector. These objects are marked with an 
"a" in the last column of Table El 



5 MEMBERSHIP 

Determining the membership status of candidates can be a 
subjective process largely dependent on the specific criteria 
set. This section details the three diagnostic tools we used 
to assign membership of the a Orionis cluster. 

5.1 Radial Velocities 

Our initial indicator of cluster membership was based on 
comparing individual RVs to the weighted mean cluster ve- 
locity. Targets lying more than 3a away from this mean were 
flagged as likely non-members (or binaries) and excluded 
from subsequent recalculations. In the 19 cases where two 
RVs were available for a given object, the average of the two 
was used in computing the ensemble velocity. The result of 
this iterative process was that 66 objects were designated as 
a Orionis RV cluster members, lying within 3a of a cluster 
mean velocity of 31.2 ± 0.2kms~^ (which has an additional 
external error of up to lkms~^). 

Our mean RV for the cluster agrees reasonably well 
with the systemic velocity found for a Ori by Morrell & 
Levato (1991) of 27kms"^; the value of 29.1 kms"^ listed 
for a Ori by Evans (1967); the peak in the RV distribution 
of low-mass PMS candidates around a Ori at 25-30 km s~^ 
found by Walter et al. (1998) and at 29.5 km s"^ by Burn- 
ingham et al. (2004). It is somewhat lower than the mean 
of 37.3 km s~^ reported for a group of candidate low-mass 
members by Zapatero-Osorio et al. (2002). 

The iterative clipping process has the potential to ex- 
clude candidate members which are cluster binary systems 
with a small enough separation to introduce RV variations 
that move them away from the cluster mean. If other in- 
formation suggests strongly enough that candidates are in- 
deed a Ori association members, then a single discrepant 
RV measurement can be sufficient to identify a short pe- 
riod, single-lined spectroscopic binary. So, all objects failing 
the RV test must be examined closely for any convincing 
signs of membership from other indicators. 

From these considerations we have 10 RV non-members 
in our sample, with heliocentric RVs ranging from -25 to 
120kms"\ If we approximate our RV selection to a band of 
width ~ 25kms~^ centred on the cluster mean and assume 
that contaminating RV non-members are spread evenly in 
RV space then we would only expect of order 1-2 contami- 
nating non-members to be left among our RV-selected mem- 
bers. 

5.2 Sodium 

Each of our a Ori candidates exhibits the gravity sensitive 
8183, 8I95A sodium doublet in their near-infrared spectra. 
The strength of this doublet is gravity sensitive (see Schi- 
avon et al. 1997). It is 2-3 times stronger in cool M dwarfs 
(M4 and later) than in PMS stars with age < lOMyr and 
similar effective temperature (e.g. Kirkpatrick et al. 1991; 
Martin, Rebolo & Zapatero-Osorio 1996; Martin, Delfosse 



& Guieu 2004). As demonstrated in Figure [7| this test may 
only be conclusive for objects with I — J > 1.6. Warmer 
objects in the a Ori cluster will have gravities only slightly 
lower than field objects with the same colour. We define an 
arbitrary empirical boundary, plotted on Figure |7| with a 
dashed line, below which we are confident that a low Nal 
EW indicates a gravity consistent with PMS status. 64 can- 
didates lie below this line and 12 above. 

The Nal EW is also very weak in giant stars (^ lA 
according to Schiavon et al. 1997), but contamination of 
our sample by M-type giants is highly improbable. An M- 
giant at the appropriate location on the CMD would have 
an absolute / magnitude < — 2 and hence would have to be 
at distances > 25 kpc. As the galactic latitude and longitude 
of a Ori are / = 206.8°, b = —17.3, such objects would be 
more than 7.4 kpc below the galactic disc, where their spatial 
density is known to be essentially zero (Branham 2003). 

5.3 Lithium 

The presence of undepleted lithium should be an excellent 
indicator of stellar youth, because it is burned in the fully 
convective interiors of low-mass PMS stars once their cores 
reach a threshold temperature. We would expect Li to start 
being depleted by more than a factor of 100 among contam- 
inating field VLMS at ages between 20Myr (at 0.35 M©) 
and 120 Myr. Field BDs with masses below about 0.06 M© 
never deplete Li, but these would have colours redder than 
I — J > 2.8 (the reddest in our sample) for ages > 140 Myr 
(Baraffe et al. 1998). According to tables of EW versus 
Li abundance presented by Zapatero-Osorio et al. (2002), 
which have appropriate temperatures for the sample con- 
sidered here (at 5 Myr 1.4 < I — J < 2.8 corresponds to 
3500 > Teff > 2650 K - Baraffe et al. 1998), undepleted 
Li should result in a Li i 6708A fine EW of about 0.5-0.6A 
falling to 0.2-0.3A after a factor of 100 depletion and becom- 
ing much weaker shortly after that. Thus the detection of 
Li I 6708A consistent with an EW of O.2A or more should be 
a very strong indicator of youth in the hottest stars of our 
sample (implying ages < 20 Myr), but is less constraining 
(ages < 140 Myr) for the coolest stars in our sample. Only 
a small fraction (~ 0.1 — 1 per cent) of any contaminating 
foreground field stars might be expected to show evidence 
for Li at this level. 

A caveat to this picture might be if there are remnants 
from an older generation of star formation in the Orion re- 
gion which might show Li. However, these older stars would 
have to lie in the foreground (to appear in the CMDs in the 
correct place) and would also have to satisfy the RV tests 
to be included as candidate members. We have to bear in 
mind that potential binary members exhibiting Li, and a 
weak Nal feature but with a discrepant RV could also be 
members of a young (10-100 Myr) foreground population. 

Figure |S] shows the 57 objects with available equivalent 
widths (and photometry) along with 19 calculated upper 
limits. We decide to rule out membership for any star which 
has EW[Li]< 0.2A at the 2-sigma level. Above this level 
the presence of Li is taken to be very strong evidence of 
cluster membership and if the Li measurement upper limit 
is inconclusive then the evidence from the other two mem- 
bership indicators is used. However, we find a number of 
objects that display weak or absent Li but which appear to 
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have low-gravity and are RV members and we will return to 
these in Section 6. 



6 DISCUSSION 

6.1 Membership Indicators 

Using three membership tests means that candidates can be 
placed into a number of categories. 

(i) Most candidates (55/76) are RV members with weak 
gravity and display evidence of significant lithium in their 
atmospheres. These are clearly the strongest cases for being 
true <J Ori members. All but one of the strong Ha emit- 
ters/accretors belong in this category. 

(ii) There are 3 candidates (targets 27, 40 and 73 - an Ha 
emitter/accretor) which have an inconclusive upper limit to 
their Li EW, but which are RV members and display low 
gravity. 

(iii) There are 3 candidates that display a low gravity, 
are RV non-members and either show Li (targets 46 and 74) 
or have an upper limit to their Li consistent with the pres- 
ence of significant Li (target 72). These are possible binary 
members of the cluster. In the case of target 72, the RV is 
demonstrably variable. 

(iv) There are 8 candidates which have no Li, are RV 
members and have a Na EW that is either inconclusive (tar- 
gets 29, 33, 34, 45, 53 and 59) or indicates a low gravity 
(targets 25 and 76). 

(v) There is 1 candidate (target 52) which shows no ev- 
idence for Li, has a discrepant RV yet a weak Na EW, in- 
dicative of low gravity. 

(vi) There are 6 candidates (targets 7, 12, 26, 28, 30 and 
38) that have very low upper limits to their Li EW, incon- 
clusive Na EW and highly discrepant RVs. These are all 
considered definite non-members 

In summary we have found what we believe are 57 firm 
members of the a Ori cluster (those displaying Li) , of which 
2 are probable short period binary systems. We have also 
identified 6 objects as clear non-members of the cluster. In- 
evitably some of our classifications are arguable, but there 
is sufficient information in Table|S|for the reader to generate 
their own classifications. 

Thirteen of our targets (an uncomfortably large num- 
ber) lie in a grey area where we find it difficult to assign 
membership, either because they show low-gravity but dis- 
crepant RVs, low-gravity but no Li, or the right RV but 
no Li. Only 2 out of these 13 show discrepant RVs, one of 
which is definitely an RV variable. Given the relatively nar- 
row range of RVs occupied by members and the spread of 
velocities occupied by the 10 (out of the total 76) targets 
with discrepant RVs, we would only have expected 1-2 non- 
members to have the "correct" RV. This presents us with 
a puzzle, because it suggests that the majority of the 11 
"grey area" objects with the "correct" RV should be cluster 
members. Three of these objects have an upper limit to the 
Li EW which could be consistent with negligibly depleted 
Li, but 8 do not, and of these, 2 (targets 25 and 76 - see 
Fig. |3(a)| l have a small Na EW that would seem to unam- 
biguously classify them as low gravity PMS objects. This 
in turn suggests there may be a problem in using a simple 
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Figure 9. /,/ — J colour-magnitude diagram. Symbols have the 
same meaning as in Fig. 1 with open squares around those objects 
found to have accretion. Evolutionary tracks for masses 0.04 — 
0.6Mq (dotted lines) and isochrones for ages 1 — 50Myr (dashed 
lines) are calculated using the models of BarafTe et al. (1998). 
The isochrones range from 1 Myr (far right) to 50 Myr and the 
evolutionary tracks are shown for masses from 0.6 Mq (at the 
top) to 0.04 M0. 



EW[Li] threshold to discriminate between cluster objects 
and non-members. On the other hand, that 2 out of 13 grey 
area objects are RV non-members contrasts with only 2 out 
of 57 Li-rich RV non-members. So unless short period bina- 
rity somehow has a tendency to reduce EW[Li] , it is probable 
that one or more of these 2 are not cluster members, despite 
appearing to be low-gravity objects (see also section 6.4); in 
which case we have the problem that EW[Na] may not be a 
foolproof membership indicator. 



6.2 A spread in Li equivalent widths? 

We have observed a large spread in the measured values of 
EW[Li]. The spread appears to be confined to those cluster 
candidates with I — J > 1.7 (see Fig. |SJ and cannot be ex- 
plained with the quoted experimental errors - many of our 
spectra are of reasonably high quality and we believe the 
statistical errors in the EW values are quite robust, if not 
pessimistic. In particular there appear to be a substantial 
fraction of the cooler stars which have EW[Li]~ O.SA. In 
addition there are also the "grey area" objects, the majority 
of which are likely to be members but most of which ap- 
pear to exhibit no Li (including three with a definite PMS 
gravity). 

It also seems that the average EW[Li] declines with in- 
creasing I — J from about 0.6A to 0.5A across our full I — J 
range. We mention this because the tables of EW versus 
Li abundance presented by Zapatero-Osorio et al. (2002) 
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show instead a slight increase in the expected EW[Li] for 
an undepleted Li abundance, of about the same magnitude. 
A plausible explanation for the decreasing average EW[Li] 
with / — J in our sample arises from the definition of the 
pseudo-continuum used to estimate EW[Li]. The continuum 
regions we have used cover a rather wide range around the 
Li feature, so they must surely contain molecular absorption 
features which become stronger at cooler temperatures. This 
would depress the pseudo-continuum leading to lower esti- 
mates of EW[Li]. 

A real Li abundance spread could be caused by deple- 
tion as Li is burned in the cores of contracting young low 
mass objects. As explained in section 5.3, factors of 100 
depletion in Li are probably needed to reduce EW[Li] to 
^ O.3A, which does not occur (according to the models of 
Baraffe et al. 1998) until ages of 20-120 Myr. Other mod- 
els, such as those of D'Antona & Mazzitelli (1997) or Siess, 
Dufour & Forestini (2000), predict similar if not slightly 
larger timescales for Li depletion. Figure |^ shows the I ver- 
sus I — J colour magnitude diagram for our targets together 
with isochrones from the Baraffe et al. (1998) models. The 
models are adjusted for the Hipparcos distance to a Ori of 
352 pc and suggest a spread of ages between < 1 and 20 Myr 
for our Li-confirmed cluster candidates. Only some of the 
oldest, most massive stars in our sample could have suf- 
fered significant Li depletion if these ages are taken literally. 
However, the "oldest" Li-poor cluster candidates among the 
bright stars in our sample also have discrepant RVs and so 
are almost certainly cluster non-members. In any case the 
apparent age spread present in the diagram is likely to be a 
vast overestimate - photometric variability, binarity and for 
some stars, the presence of discs, can perhaps entirely ex- 
plain such spreads (e.g. Ifartmann 2001). We conclude that 
Li-depletion among members of the a Ori cluster is unlikely 
to be the cause of spreads in EW[Li] unless current PMS 
evolutionary models have very serious flaws. 

An alternate possibility is that some of our cluster can- 
didates are genuinely depleted of Li because they are older 
and are not cluster members. They could belong to a fore- 
ground population of older stars that are at an appropriate 
distance to place them close to the a Ori cluster locus in 
CMDs (see Alcala et al. 2000; Walter et al. 2000). Such ob- 
jects would need to have a similar RV to the cluster (which is 
quite different to the general fleld population) , be old enough 
to have depleted a factor of 100 or more of their initial Li 
and be yet be young enough to exhibit weak Na lines. There 
may be a vanishingly small range of ages that can satisfy 
these requirements and whilst a slightly older population 
(the Orion OBla subgroup; age 11 Myr) does exist in some 
parts of the Orion star forming region, it is at a similar dis- 
tance and there is no evidence of this older component in 
the direction of a Ori (Sherry 2003). To conclude, we find 
this possibility unlikely. 

A second class of explanations are those in which the 
spread is apparent in the EWs because of deficiencies in 
the atmospheric analysis or the data reduction. The EW- 
abundance relationships we are using to interpret our results 
show that the EW[Li] versus Tcff relationship is quite fiat 
over the Tefs range considered here. It is therefore unlikely 
that uncertainties in the colours play any significant role in 
the EW[Li] scatter. Likewise, temperature inhomogeneities 
on the stellar surfaces caused by magnetic fields are are also 



unlikely to have much effect. Some photospheric lines can be 
veiled in young stars that are undergoing accretion, but the 
majority of stars in our sample are accreting at a rate M < 
lO~^^M0/year (see section 6.5) which is insufficient to cause 
such an effect (MuzeroUe et al. 2003). Also, the 5 or 6 objects 
which do show evidence for ongoing accretion have EW[Li] 
values of 0.32A to 1.19A, which are not systematically low 
compared with the rest of the sample. 

We have also considered possible data reduction prob- 
lems such as an uncertain and perhaps systematically under- 
estimated scattered light or sky contributions which would 
preferentially affect the fainter targets. The level of scat- 
tered light in the spectrograph is simply too small to be of 
concern here. The sky contribution is more significant, but 
we estimate that it is subtracted with a precision of better 
then 5 per cent based on the scatter in the fibre-to-fibre rel- 
ative efficiencies determined from three successive offset sky 
exposures. We have confirmed that this would not signifi- 
cantly increase the quoted EW[Li] uncertainties for any of 
our targets. 

Finally we note that a spread in EW[Li] among VLMS 
and BDs has been seen before, although not extensively dis- 
cussed. Both Jorgens & Guenther (2001) and Natta et al. 
(2004) have reported significant EW[Li] spreads of < 0.09 
to O.63A and 0.2 to 0.5A respectively, among a small group 
of very low mass objects in the young Chamaeleon I clus- 
ter. It is also notable that a subset of these objects were 
observed earlier by Neuhauser & Comeron (1999), albeit at 
lower resolution, who found quite different EW[Li] for some 
objects. 

We can reach no firm conclusion here. The EW[Li] 
spread demands further attention with better spectra, if only 
to rule out the possibility of major problems with the cur- 
rent evolutionary models. Our results cast doubt on any 
selection technique which rules out membership of young 
clusters solely on the basis of a small EW[Li]. 

6.3 Contamination and the mass function 

Previous work on the a Ori cluster has attempted to esti- 
mate the mass function of VLMS and BDs using photomet- 
rically selected samples, finding dn/dlogm oc m~°', with 
Q = -0.2 ± 0.4 for 0.01 < m < 0.2 Mq (Bejar et al. 1999, 
2001). Only a minority of this sample have been followed 
up spectroscopically and in most cases, this spectroscopy is 
limited to obtaining a spectral type and noting the presence 
of moderate to strong Ha emission in low resolution (20A 
FWHM) spectra. We contend that such an approach may 
not yield reliable results. As we discussed earlier, the main 
source of any contamination in the sample will be late-type 
foreground M dwarfs that have spectral types, colours, and 
in many cases Ha emission, that are similar to members of 
(J Ori. A very large Ha EW of perhaps 20-40A or more is 
associated with accretion and hence youth, but is present 
in only a small fraction of a Ori members (see section . 
Measurement of radial velocities, Li and weak gravity using 
Na and K line strengths will be much more reliable, but is 
only possible with higher resolution spectra. 

Arguments based upon the space density of field M 
dwarfs, have suggested that any contamination in a pho- 
tometrically selected sample may be quite low (see Bejar 
et al. 1999 and below). Our spectroscopy and membership 
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classification partially reinforce this view. We have found at 
least 57 (and probably closer to 64) likely members from a 
sample of 76 selected only on the basis of their / magnitudes 
and R ~ I colours. Cluster members are found close to the 
edge of the sample distribution in both the I ys R — I and 
7 vs / — J CMDs, suggesting that the photometric selection 
criteria may have been too restrictive to include all the pos- 
sible members of a Ori. RV measurements of another faint 
sample of a Ori candidates spread over a wider region of 
the 7 vs i? - 7 CMD by Burningham et al. (2004) suggest 
that few members lie outside the region we have considered. 
However, because of this uncertainty and because we have 
not observed a complete sample of photometric candidates 
over our restricted field of view, we limit ourselves in this 
paper to asking whether a mass function derived solely by 
photometric selection is likely to be significantly affected by 
non-member contamination? 

The definite non-members we have found are at the 
brighter end of the sample (14.8 < 7 < 16.3) - where the 
contamination level is 15-25 per cent (depending on the sta- 
tus of the "grey area objects"). Based solely on their low 
EW[Na], it is possible that all of the fainter objects are 
members, but the contamination level could be as high as 
20 per cent if all the "grey area" objects turn out to be 
non-members. Contamination at this level would have min- 
imal impact on the value of a, especially when compared to 
the effect of statistical, distance and model-dependent un- 
certainties. 

To check our estimates of contamination we have simu- 
lated an 7 versus R — I diagram populated by field dwarfs. 
The details of this simulation, which uses the local field 
dwarf luminosity function and space density, are presented 
in Jeffries et al. (2004). We simulated a 0.5 square degree 
field towards o Ori, similar to that area from which our spec- 
troscopy targets were chosen. An I versus 7? — 7 diagram was 
generated for objects up to 500 pc distant. We then selected 
objects with 14.8 < 7 < 18.2 and in a strip 0.6 magnitudes 
wide in 7? — 7 and centred on the 5 Myr isochrone shown in 
Fig 1. We found that 41 simulated M-dwarfs lay within this 
strip. In the real photometry dataset there are 115 objects in 
the same strip, of which 74 were observed spectroscopically 
and 57 — 70 are found to be members depending on the sta- 
tus of "grey area" objects. Hence we would expect to have 
observed about 25 non-members in our spectroscopic sam- 
ple, compared to the 4-17 we have found. We consider this 
reasonable agreement, although any discrepancy might be 
explained by a tendency for us to prefer to allocate fibres to 
spectroscopic targets which lay close to the 5 Myr isochrone. 
Most of the contaminants will lie on the blueward side of the 
selected strip. Indeed only 2 simulated contaminants lay on 
the red side of the 5 Myr isochrone and Fig. 1 shows that 
all the confirmed non-members and all but one of the "grey 
area" objects lie blueward of the 5 Myr isochrone. 

It is interesting to note that had we applied a strict 
age^ 20 Myr criterion using the 7 vs 7 — J diagram shown 
in Fig. |5] we would have excluded all of the definite non- 
members from our sample. The 4 "grey area" objects that 
would be excluded on this basis have no Li, have an am- 
biguous EW[Na] and in two cases an RV only marginally 
consistent with cluster membership. Selection on the 7 ver- 
sus I — J diagram seems therefore to have the potential for 
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Figure 10. The efficiency with which we are able to detect bi- 
naries as a function of separation. The solid curve is averaged 
over the 57 objects in sample A (see text). A similar curve for 
sample B is almost indistinguishable. We have assumed circular 
orbits, random inclination angles, a flat mass ratio distribution 
for 0.2 < g < 1 and primary masses taken from Fig. 1^ 

being more discriminating in choosing a sample of genuine 
cluster candidates. 

6.4 Binarity 

The data we have obtained offer us the opportunity to inves- 
tigate the fraction of VLMS and BDs in the a Ori cluster 
which are in short period binary systems. A test can be 
formulated to see whether the RVs of cluster members are 
(a) constant, if measured more than once, and (b) consis- 
tent with the mean cluster RV. The fraction of short period 
binaries identified in this way can then be compared with 
models of the binary fraction as a function of separation, 
a, and mass ratio, g, that also take into account the RV 
sampling and uncertainties. 

A difficulty here is deciding which candidate VLMS 
and BDs to include as genuine members. Clearly we can- 
not use measured RV as a criterion! We make two choices, 
one more restrictive than the other. The first sample is all 
those VLMS/BD candidates which show clear evidence for 
Li in their spectra - a total of 57 objects, which we will call 
sample A. The second sample, includes the first as a subset, 
but also includes all those objects with an EW[Na] which 
would seem to unambiguously classify them as low-gravity 
objects - a total of 64 objects, which we will label as sam- 
ple B. Since these bracket a range of possibilities for the 
binary fraction, we perform calculations using both. 

We applied a test to the radial velocities of the can- 
didate cluster members in each sample in order to identify 
binary systems. In effect, we were testing the goodness-of-fit 
of the mean cluster RV to the observed RV(s), adding the 
uncertainty in the cluster mean in quadrature to each error 
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Figure 11. On the left we show the assumed distribution of binary separations. The top plot has a peak at 15 au and a Gaussian sigma 
of 1.53 in logarithmic units. They have been normalised to show the approximate binary fractions seen in sample of nearby M-dwarfs 
(see Fischer &; Marcy 1992, figure 2), with a total binary fraction of 0.42. The bottom panel shows a flat distribution with a similar 
normalisation. On the right we show the corresponding probability that these distributions produce the observed numbers of binary 
systems. The solid curves show the results for sample A (2 binaries from 57), and the dashed curves show the results for sample B (4 
binaries from 64). 



bar. The threshold was set so that there was a probabil- 
ity of < 10~^ of a constant RV cluster member failing the 
test. Two objects from sample A (targets 46 and 74) were 
identified as binary and a further 2 objects (targets 52, and 
72) in sample B. Then, assuming circular orbits and primary 
masses taken from Fig. |^ we used the mean times of obser- 
vation and the RV error bars for each target to calculate 
the probability that a binary RV variation would have been 
detected at this level, both as a function of a and q and 
averaged over all inclination angles. Summing these proba- 
bilities over all targets and over particular distributions of q 
yields a "detection efficiency" as a function of binary sepa- 
ration for that particular sample. We can then calculate the 
probability that a given binary fraction, q and a distribution 
can result in the number of binaries observed (see Maxted 
et al. 2001). 

As an example. Fig. 1101 shows how efficiently (on aver- 
age) we can detect binary systems for a given separation for 
samples A and B. In the absence of any definitive informa- 
tion on the distribution of q in VLMS /BDs we assume a flat 
meiss ratio distribution for 0.2 < q ^ 1 which is reasonably 
consistent with what is observed for M-dwarfs (Fischer & 
Marcy 1992). If the q distribution were peaked towards q = 1 



then binaries become more detectable and the required bi- 
nary fraction to explain a set of observations becomes lower 
(and vice- versa). It is clear from this plot that our survey 
is not sensitive to systems with separations of log a > (in 
au) . Conversely our detection efficiency for closely separated 
binaries is very high. 

Given the sparse nature of our RV dataset we limit our- 
selves to an attempt to answer one simple question. Is the 
number of candidate binary systems we have observed con- 
sistent with the overall binary fraction (42 ± 9 per cent) and 
Gaussian log a distribution proposed by Fischer & Marcy 
(1992) for M-dwarfs within 20 pc? In one sense we already 
know that for BDs the answer to this is no. There are too few 
field and cluster BD-BD binary systems with a > 20 au com- 
pared with M-dwarf binaries in the same separation range, 
(e.g. Martm et al. 2003; Gizis et al. 2003). However, we are 
probing the small separation end of this distribution which 
is a separate problem. We begin by using Fig. 1101 to calcu- 
late the probability that a given binary fraction would have 
resulted in the number of observed binary detections, as a 
function of log a. We then convolve this with an assumed 
separation distribution that is a Gaussian with a peak at 
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log a = 1.16 and a sigma of 1.53 (in au)^. The normalisation 
of the function is effectively set by the number of binaries we 
have observed. A plot of probability versus binary fraction 
can then be produced by integrating over any given range 
of log a, although as the detection efficiency becomes neg- 
ligible for log a > (see Fig. Il()|l we confine ourselves to 
-3 < loga < 0. 

In the case of sample A and assuming that the Gaussian 
log a distribution shown in Fig. Illb is appropriate, we find 
that the most likely binary fraction (between —3 < log a < 
only) is 0.19, with a lower limit of > 0.13 at a 90 per 
cent confidence level (see Fig. lllb '). This is inconsistent with 
the estimate of about 5-7 per cent binarity for this range 
of log a (Marcy & Benitz 1989; Fischer & Marcy 1992) in 
nearby M-dwarfs. For sample B where there is an even higher 
observed binary fraction the probability curve peaks at a 
binary fraction of 0.47 and with a lower limit of > 0.17 at a 
90 per cent confidence level. 

These results imply that if VLMS /BDs do have a sep- 
aration distribution that decreases towards smaller separa- 
tions in a similar way to field dwarfs, then we have observed 
too many to be consistent with the small binary fraction 
found for field M-dwarfs at small separations. This in turn 
means that either the distribution of separations is quite 
different in the VLMS/BDs and/or the binary fraction at 
small separations is very high. To test the sensitivity to the 
assumed distribution of log a, we have performed another 
calculation assuming a flat distribution (see Fig. Illb ). In 
this case sample A returns a most probable binary fraction 
(-3 < loga < 0) of only 0.07, with a lower limit of 0.07 at 
90 per cent confidence. So, depending on the detail of the 
separation distribution it seems that fraction of VLMS /BDs 
in short period binary systems could be quite small (if the 
log a distribution is much flatter than found for higher mass 
stars) or could certainly be as high as the 0.5 suggested by 
Pinfield et al. (2003) for VLMS/BDs in the Pleiades. Ei- 
ther way the binary properties for short-period VLMS/BDs 
would be quite different from those of M-dwarfs, as indeed 
they are at larger separations. 

There are several caveats to this interesting conclusion, 
(i) It is based on very small number statistics, which al- 
though statistically taken account of in our calculations, 
render the result vulnerable to contamination of the binary 
sample, (ii) Greater consistency with the M-dwarf distri- 
bution would be achieved if some of the candidate binary 
members were not binaries. This might be the case if our 
RV errors were underestimated, but in that case we would 
be less sensitive to binaries and the deduced binary frac- 
tions would remain similar, (iii) The binary fraction would 
be smaller if some of the binary candidates turn out not 
members of the cluster. This certainly applies to sample B 
which contains two "grey area" binary candidates for which 
we are not certain of cluster membership. Indeed the large 
implied binary fraction could be viewed as evidence that this 
is the case. However our main conclusion holds even if we 



^ Fischer & Marcy favoured a similar period distribution to that 
of G-dwarfs found by Duquennoy & Mayor (1991). Here we have 
done the same which implies that the peak position of the a distri- 
bution should be reduced by a factor of two, from 30 au to 15 au, 
to account for the lower system mass among VLMS/BDs 



consider just the sample in which Li was detected, (iv) If 
the q distribution for VLMS/BDs was sharply peaked to- 
wards q — 1 this would make it easier to detect binaries and 
thus our efficiency would be greater and the required binary 
fraction lower. 

Our results seem reasonably consistent with the survey 
of 26 field VLMS/BDs conducted by Guenther & Wuchterl 
(2003). They used 2-3 RV measurements with precision 0.1- 
lkms~^ spread over 1-2 months. They find 3 binaries in 
this sample and hence deduce a binary fraction of at least 
12 per cent for periods of < 100 days (roughly equivalent to 
0.2 au). If we perform our calculation for —3 < loga < —0.7 
for the case of sample A and the Fischer & Marcy (1992) 
distribution, we obtain a most probable binary fraction of 11 
per cent. Both of these numbers should be compared with 
the ~ 3 per cent binarity implied by the Fischer & Marcy 
analysis over a similar separation range. To make further 
and more decisive progress would require a multi-epoch RV 
survey of a large sample of VLMS /BDs in the a Ori cluster 
and elsewhere. Observations spread over several months and 
with a factor of a few better precision would give sensitiv- 
ity to separations as large as the minimum observable with 
adaptive optics systems for nearby BDs and would allow 
some estimation of the log a distribution. 

6.5 Disc Fractions 

The frequency and lifetimes of discs around young, 
VLMS/BDs may be an important diagnostic tool in identify- 
ing the dominant (sub)stellar formation scenario. A shorter 
disc lifetime for lower mass objects as a result of dynamical 
ejection of stellar embryos from their initial gas reservoirs 
is qualitatively supported by numerical models (Bate et al. 
2003). Haisch et al. (2001) determined disc fractions (using 
K — L near infrared excess) for low-mass stars (0.2-1 M©) 
in several clusters, demonstrating a rapid decrease in disc 
frequency with age - initial frequencies may be as high as 
~ 80 per cent, with approximately half of the discs being lost 
within < 3 Myr, and an overall lifetime of ~ 6 Myr. Thus, at 
an age of 3 — 7 Myr some dissipation of circumstellar discs 
is certainly expected within the a Ori cluster. Using simi- 
lar diagnostics, Oliveira, Jeffries & van Loon (2004) found 
that ~ 50 per cent of low-mass stars around a Ori possessed 
circumstellar material, broadly agreeing with this scenario. 
Liu et al. (2003) and Jayawardhana et al. (2003b) present 
K — L measurements for BDs in a number of clusters which 
support a similar timescale for disc dissipation in lower mass 
objects. 

Five or possibly six of our very low-mass o Ori clus- 
ter members show evidence for the presence of accretion 
on the basis of broad Ha emission. Three of these objects 
are probably substellar according to Fig. |5| The deduced 
accretion frequency of 10 ± 5 per cent is much lower than 
the frequency of discs found among higher mass objects by 
Oliveira et al. (2004). However, different disc indicators can- 
not be considered equivalent. Natta et al. (2004) show that 
a broadened Ha profile with a full width at 10 per cent 
of maximum of > 270kms~^ is sensitive to accretion rates 
> 10"" Mq yr"^ (for VLMS/BDs), whereas a if- L excess 
can be seen even from a passive (reprocessed light only) disc 
with M > IQ-'^ Mq (Wood et al. 2002; Walker et al. 2004). 
For example, 23 of the brighter objects in our sample (all 
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with / < 16.75) have K — L measurements in OUveira et al. 
(2004) or Ohveira et al. (in preparation). Of these, 6 to 9 
(depending on the exact criteria) show a. K — L excess, but 
only targets 10, 13 and 37 have been identified as accretors. 
Jayawardhana et al. (2003b) also find that 2 out of 6 BDs in 
0" Ori have K — L excesses, but we have found no evidence 
for accretion in either of these (S Ori 12 and S Ori 40). 

A more sensible comparison would be with the fre- 
quency of strong, accretion-related Ha emission in stars of 
higher mass in the a Ori cluster and VLMS /BDs in younger 
and older clusters. Among younger clusters, Jayawardhana 
et al. (2003a) find the fraction of accreting substellar objects 
(based on the width of the Ha line) to be about 50 per cent 
for ages 2 Myr (in IC348). MuzeroUe et al. (2003) mea- 
sured a small sample of eight a Ori candidates, finding no 
accretors. S Ori 12, 29 and 40 are in common with this study 
and we also find no evidence for accretion in these objects. 
Neither do they find any VLMS/BD accretors in the older 
(~ 8 Myr) Sco-Cen association. Our measurements from a 
much larger sample clarify this trend and pin down the ac- 
cretion timescale for most VLMS/BDs to be less than the 
age of the a Ori cluster. 

The Hq results for VLMS/BDs in the a Ori cluster are 
similar to those for higher mass stars. Wolk (1996) states 
that 12 per cent (6 out of 49) among a sample of o Ori 
candidates with spectral types KO-MO are classical T-Tauri 
stars with accretion discs, using EW[Hq] > lOA as his cri- 
terion. Using the same Ha criterion, Zapatero-Osorio et al. 
(2002) find one star from 11 shows accretion in the same 
spectral type range. Among younger clusters, the prepon- 
derance of accretion-induced Ha emission is much higher. 
For example, 63 per cent (30/48) stars with M > 0.2 Mq 
are accretors in the Taurus- Auriga (age 1-2 Myr) sample of 
White & Ghez (2001). 

Hence, it seems that, similarly to the K ~L disc indica- 
tor, the fraction of objects with discs betrayed by their ac- 
cretion signatures declines at a rate that is roughly indepen- 
dent of mass, but that the decline takes place more rapidly. 
Interpretation of this result is complicated by the mass- 
dependence of the sensitivity of these tests. Most low-mass 
classical T-Tauri stars accrete at rates that are orders of 
magnitude higher than the threshold below which their ac- 
cretion would be undetectable. Accreting BDs on the other 
hand are detected at rates right down to a smaller thresh- 
old. This suggests that the transition between accreting and 
non-accreting status may happen quite rapidly for low-mass 
stars and that the fraction of accreting BDs may be an un- 
derestimate. Either way it seems that accretion may become 
undetectable even whilst a significant mass of disc material 
is still present. If accretion at a rate ~ 10~^^ Mq yr~^ is 
present in BDs at an age of ~ 5 Myr, it implies disc masses of 
order 5 x 10~^ Mq. The recent models of flared discs around 
BDs by Walker et al. (2004), do indeed suggest that near-IR 
excesses should be present for these and even less massive 
discs. 

The presence of accretion discs in a significant fraction 
of VLMS/BDs at an age of a few Myr and in similar num- 
bers to those seen in higher mass stars at the same age, 
does not support the idea of VLMS/BDs forming by a dy- 
namical ejection process in which their discs are truncated. 
Instead it argues for a common formation mechanism for 
low-mass stars, VLMS and BDs. However, whether such ob- 



servations can definitively rule out the ejection model must 
await detailed predictions of the initial BD disc properties 
and masses from such a scenario (e.g. see the discussion in 
Liu et al. 2003 and Natta et al. 2004). 



7 CONCLUSIONS 

We have used intermediate resolution spectroscopy to refine 
a set of candidate VLMS and BD members of the a Ori 
cluster (age ~ 3 — 7 Myr) that were chosen on the basis of 
their RI photometry. We have measured radial velocities, 
the strengths of absorption lines due to neutral lithium and 
sodium and searched for signs of accretion activity. 

We find that the majority of the photometrically se- 
lected candidates are genuine cluster members on the basis 
that they exhibit an Li I 6708A absorption feature, have ra- 
dial velocities consistent with cluster membership and have 
a Na I 8183/8195A absorption doublet that is weaker than 
main sequence M-dwarfs of a similar spectral type. The clus- 
ter members we have identified have model- and distance- 
dependent masses of 0.03 < M < 0.35 Mq. Approximately 
23 of our targets are likely to have substellar masses when 
measured in this way. 

We have found another class of candidate where the ev- 
idence is ambiguous. These objects either have upper limits 
to their Li line EWs that may be consistent with cluster 
membership, have radial velocities that indicate they are 
non-members or cluster binaries or have no Li but radial ve- 
locities consistent with cluster membership. Several of these 
objects have weak Na I absorption that indicates a youthful 
status. These objects, together with the apparent spread in 
the strength of the Li I 6708A absorption EW among cluster 
members, lead us to conclude that there may be deficiencies 
in our understanding of the formation of this feature in cool, 
low gravity objects. It may be unsafe to rule out the mem- 
bership of a VLMS or BD from a young cluster or association 
solely on the basis of a weak or absent Li feature. 

Despite these problems we find that overall, photomet- 
ric selection alone is reasonably effective in isolating mem- 
bers of the a Ori cluster. We estimate contamination rates 
of 15-25 per cent among the higher mass stars of our sample 
and < 20 percent at the lower mass end. Contamination at 
this level will not significantly change the form of the cluster 
mass function that has been deduced from photometrically 
selected samples. However, we do show that selection in the 
I versus I — J diagram is potentially more discriminating 
than the / versus R — I diagram. 

We have identified 4 candidate binary cluster members 
on the basis of a varying RV or an RV discrepant from the 
cluster mean combined with other indications of member- 
ship. Two of these candidates show signs of Li in their spec- 
tra. A careful consideration of the RV errors and sampling 
in our data enables us to determine how efficiently we can 
detect binary systems as a function of primary mass, mass 
ratio and separation. Combining this calculation with the 
number of binaries we have observed and plausible distribu- 
tions of binary mass ratio and separation, we determine that 
we have found too many binary systems to be accounted 
for by binary properties similar to those in higher mass 
stars. Instead, the short-period binary frequency could be 
much larger in our VLMS/BDs (perhaps as high 0.5 for 
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—3 < log a < 0) or the binary frequency is lower but the 
distribution with log a must be flatter than seen in higher 
mass field M-dwarfs. Caveats to this conclusion must be con- 
sidered. The results are crucially dependent on the correct 
identification of a small number of binary members of the 
cluster. However, these conclusions still apply even if we re- 
strict ourselves to considering cluster members where the 
presence of Li has been positively identified. Further work 
will be needed to hone this result - a multi-epoch RV survey 
of a large sample of a Ori VLMS/BDs with a better RV pre- 
cision has the capability to ascertain the binary frequency 
and separation distribution at separations smaller than ac- 
cessible by adaptive optics observations of nearby BDs. 

A search for broadened Ha emission features associ- 
ated with a circum (sub) stellar accretion disc has revealed 
5 or 6 accreting objects, 3 of which are found among the 
substellar candidates. The fraction of VLMS/BD accretors 
(10 ±5 per cent) is similar to that found among higher mass 
stars in the cluster using similar diagnostics, but much lower 
than the fraction found amongst low-mass stars, VLMS and 
BDs in younger clusters. The timescale for accretion rates 
to drop to < 10~^^ M0yr~^ is thus constrained to be less 
than the age of the a Ori cluster for most VLMS/BDs. The 
similar timescale for a decline in accretion rates for both 
VLMS/BDs and higher mass stars does not support the idea 
that VLMS/BDs have truncated, lower mass discs as a result 
of dynamical ejection. However, in the absence of detailed 
predictions of disc masses and properties for this scenario, it 
cannot be completely ruled out. Instead the observations are 
consistent with the notion that low-mass stars and BDs form 
a continuum of outcomes from a single formation process in 
which the presence of circumstellar material is a natural 
consequence. 
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Table 6. Positions, I magnitudes, I — J and R — I colours, heliocentric radial velocities (nights 1 and 2) and measured equivalent widths 
for our targets. Photometric uncertainties in / and R — I are about ±0.04 and ±0.05 respectively (see section 2). J magnitudes are from 
2MASS. Alternate names from Bejar et al. (2001) are given. The final three columns indicate our classification of membership according 
to whether the EW[Na] indicates a low gravity (Y) or is ambiguous (?); whether the RV is consistent with membership and finally 
whether an object is classed as a member (Y), non-member (N) or is in a grey-area (?). Objects with broad Ha emission or which are 
cluster binary candidates are labelled with an a (for accretor) or b. 
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